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ABSTRACT 

Context. Blazars are astrophysical sources whose emission is dominated by non-thermal processes, typically interpreted as syn¬ 
chrotron and inverse Compton emission. Although the general picture is rather robust and consistent with observations, many aspects 
are still unexplored. 

Aims. Polarimetric monitoring can offer a wealth of information about the physical processes in blazars. Models with largely different 
physical ingredients can often provide almost indistinguishable predictions for the total flux, but usually are characterized by markedly 
different polarization properties. We explore, with a pilot study, the possibility to derive structural information about the emitting 
regions of blazars by means of a joint analysis of rapid variability of the total and polarized flux at optical wavelengths. 

Methods. Short timescale (from tens of seconds to a couple of minutes) optical linear polarimetry and photometry for two blazars, 
BL Lacertae and PKS 1424+240, was carried out with the PAOLO polarimeter at the 3.6 m Telescopio Nazionale Galileo. Several 
hours of almost continuous observations were obtained for both sources. 

Results. Our intense monitoring allowed us to draw strongly different scenarios for BL Lacertae and PKS 1424+240, with the former 
characterized by intense variability on time-scales from hours to a few minutes and the latter practically constant in total flux. Es¬ 
sentially the same behavior is observed for the polarized flux and the position angle. The variability time-scales turned out to be as 
short as a few minutes, although involving only a few percent variation of the flux. The polarization variability time-scale is generally 
consistent with the total flux variability. Total and polarized flux appear to be essentially uncorrelated. However, even during our 
relatively short monitoring, different regimes can be singled out. 

Conclusions. No simple scenario is able to satisfactorily model the very rich phenomenology exhibited in our data. Detailed numerical 
simulations show that the emitting region should be characterized by some symmetry, and the inclusion of turbulence for the magnetic 
field may constitute the missing ingredient for a more complete interpretation of the data. 

Key words. Polarization - BL Lacertae objects: BL Lacertae - BL Lacertae objects: PKS 1424+240 


1. Introduction 

Blazars, the subclass of active galactic nuclei (AGN) showing 
jets almost aligned with the observer’s line of sight (Blandford 
|&Rees|1978||Urry &Padovani|1995| l, offer an invaluable labora¬ 


* Partly based on data obtained at the INAF / Telescopio Nazionale 
Galileo at the Canary Island of La Palma under program Id: 
A29TAC_21 (PI: S. Covino). 


tory of physics. Their spectral energy distribution (SED) shows 
a characteristic double-hump shape and is usually well modeled 


as due to synchrotron and inverse Compton radiation (Ghisellini 


|et al.|1998) . Relativistic Doppler boosting of the observed emis 
sion is likely involved in the large-amplitude variability observed 
essentially at all frequencies (e.g. Ghisellini et al.|l993jl. 


Although the interpretative scenario seems to be well estab¬ 
lished, many open problems are still present. The availability 
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of continuously improving multi-wavelength (MW) data has re¬ 
vealed the need for more sophisticated approaches, with models 
assuming that the observed emission originates in multiple zones 
with typically independent physical parameters (e.g. Aleksic et 
[aT7||2Q12] >. The possibility of inhomogeneity in the emitting re¬ 
gion, mimicked by multi-zone models, is definitely plausible. 
However, this immediately introduces a strong degeneracy in the 
already large parameter space, in turn requiring additional in¬ 
formation to disentangle the various possible components in the 
observed emission. 

The dominance of non-thermal emission processes (e.g. syn¬ 
chrotron radiation, etc.) in the blazar emission suggests that a 
wealth of information might come from polarimetric studies 
([Larionov et al.||2013| Sorcia et al. [2013. Sasa da et al.||2014] 


|Sorcia et~ 2014 Zhang et al.|2014 Itohetal. 2015 to mention 

some of the most recent papers). In the optical, the detection 
of polarized emission was considered the smoking-gun signa¬ 
ture for synchrotron emission from a non-thermal distribution of 
electrons) Angel & Stockman] 1980) 1. In general, the addition of 
polarimetric data to the modeling of blazar photometric/spectral 
information has widely shown its potential to derive information 
about, e.g., the magnetic field state (e.g. Lyutikov et al. 2005; 


Marscher 


2014|l, or to drive the modeling of different SED com 


ponents (Barres de Almeida et al. 2014 |, 


A relatively less explored regime is that of short timescale 
polarimetry (|Tommasi et ah 2001 a|b[ Andruchow et al.1[2005; 


|Sasada et al.|2008||Chandra et al.|26l2[|Itoh et al ,|20 1 3[ >. Short 

timescale photometry, on the contrary, is indeed a common prac¬ 
tice in the field and has revealed to be a powerful diagnostic 


technique (Montagni et al. 

2006;|Rani et al. 2010} |Danforth et 

al. 20131 Zhang et al. 2013 

Sandrinelli et al. 2014), also in the 

very-high energy regime (e.g. Aharonian et al. 2007; Albert et 

al. 20071 Abdo et al.|2010||Foschini et al.|2013 ). 


In this paper we present and discuss well-sampled obser¬ 
vations of two blazars: BLLacertae (hereinafter BLLac) and 
PKS 1424+240. The observations were carried out with the opti¬ 
cal polarimeter PAOL( 0 equipping the 3.6 m INAF / Telescopio 
Nazionale Galileo (TNG) at the Canary Island of La Palma. The 
relatively large collective area of the TNG enabled us to explore 
time scales as short as several tens of seconds in both photometry 
and polarimetry. 

The paper is organized as follows: observations are described 
in Sect.[2] In Sect.[3]results of the analyses and a general discus¬ 
sion are presented, and conclusions are drawn in Sect. [4] 


2. Observations 

PAOLO is an optical polarimeter integrated in the Naysmith fo¬ 
cus instmment DOLORE50 at the TNG. The observations pre¬ 
sented here were part of the commissioning and scientific activ¬ 
ities of the instmment. 

BLLac is the prototype of the class of BLLac objects, and is 
located at a redshift z = 0.069 ( [Miller & Hawley \99l\. The host 
is a fairly bright and massive elliptical galaxy (|Scarpa et aT 2000 


|Hyvonen et al.|2007) . Due to its relative proximity it is one of the 
most widely studied objects of the class. PKS 1424+240 is also 


a BLLac object and its redshift is still uncertain. Furniss et al. 
( |2013| > report a lower limit at z > 0.6, which can make it one of 
most luminous objects in its class. Its host galaxy was possibly 
detected by Meisner & Romani (2010) at typically a few percent 


of the nuclear emission, although jScarpa et al.| ( ]2000j i reported 
much fainter limits. 

BL Lac was observed for about 8 hours during the night of 
2012 September 1-2. The observations consisted of short in¬ 
tegrations of about 20-40 s each with the r filter, interrupted ev¬ 
ery ~ 45 min to observe polarized and unpolarized polarimetric 
standard stars (BD+28d4211, W2149+021, HD 204827) for a 
total of more than 300 data points. The data reduction is carried 
out following standard procedures and aperture photometry is 
performed using custom too 1:0 Photometric calibration was se¬ 
cured by comparison with isolated unsaturated stars in the field 
with magnitudes derived by the APASS catalogue^] Photometric 
and polarimetric light curves are shown in Fig.[T| 

PKS 1424+240 was observed for about 5 hours during the 
night of 2014 June 1-2. The observations consisted of short 
integrations of 1-2 min each with the r filter interrupted at the 
beginning and at the end of the sequence to observe an unpolar¬ 
ized polarimetric standard star (GD 319) for a total of more than 
100 data points. Reduction and calibration were carried out as 
for BLLac. Photometric and polarimetric light curves are shown 
in Fig. [2] 

The removal of the few percent instrumental polarization 
typical of Nasmyth focus instruments (Tinbergen 2007; Witzel et 
al.|201 1| Covino et al.|2014] i can be carried out rather efficiently 
and with PAOLO we can estimate ( [Covino et al.|2014| a resid¬ 
ual r.m.s. of the order of ~ 0.2% or better. If the observations 
cover a limited range in hour angles the correction is generally 
more accurate. This is a systematic uncertainty superposed onto 
our observations and it is already included in the reported errors 
for our data. The results reported here supersede the preliminary 
ones shown in Covino et al. |( |2014] l. 

Where required, x l minimization is performed by using 
the downhill (Nelder-Mead) simplex algorithm as coded in the 
pythor0 scipy. optimiz01ibrary, v. 0.14.0. The error search 
is carried out following |Cash| ( |T976) l. Throughout this paper the 
reported uncertainties are at lcr. 

Distances are computed assuming a ACDM-universe with 
Da = 0.73, Q m = 0.27, and Ho = 71 kms -1 Mpc~' (Komatsu et 


|aT7||2011| ). Magnitudes are in the AB system. Flux densities are 
computed following Fukugita et al. ( 1996) . The raw and reduced 
data discussed here are available from the authors upon request. 


3. Results and discussion 

BLLac and PKS 1424+240 are sources belonging to the same 
class and, during our observations, also showed a comparable 
brightness. This is already a remarkable finding since the latter is 
more than one order of magnitude farther away than the former. 
PKS 1424+240 is therefore intrinsically about 100 times more 
luminous in the optical than BLLac in the considered period. 
The host galaxy of BLLac was measured at R ~ 15.5 (Scarpa] 
[et al.||2000[ ), roughly 30% of the source luminosity during our 
observations. The source showed intense short-term variability, 
as expected for a blazar which has previously been found to be 
strongly variable at any time-scale ( [Raiteri et al.|2013| l. On the 
contrary, PKS 1424+240 was remarkably stable during the ob¬ 
servations with slow (hours) variations at most at a few percent 
level. This behavior is rather unexpected although this source 
presented a less intense variability (at least compared to BLLac) 


https://pypi.python.org/pypi/SRPAstro.FITS/ 


1 http://www.tng.iac.es/instruments/lrs/paolo.html 

2 http://www.tng.iac.es/instruments/lrs/ 


http://www.aavso.org/apass 

5 http://www.python.org 

6 


http://www.scipy.org/SciPyPackages/Optimize 
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Fig. 1. PAOLO observations of BLLac. In the top panel we show the 
magnitude of the source (AB mag) not corrected for Galactic reddening 
and for the host galaxy brightness. In the middle panel we show the 
polarization degree and in the bottom panel the position angle. 


Fig. 2. PAOLO observations of PKS 1424+240. The top panel shows 
the magnitude of the source (AB mag) versus time not corrected for 
Galactic reddening and for the host galaxy brightness. The middle panel 
shows the polarization degree and the bottom panel the position angle 
versus time. 


during long-term monitoring campaigns (e.g. Archambault et al. 
2014; Aleksic et al.|2014| > and in particular close to our observa¬ 
tion epodj^J 


3.1. Analysis of flux variability 

The rapid variability observed in BLLac, although in most cases 
of rather low level in absolute terms (~ 5 - 10%), is character¬ 
ized by a fair number of well sampled rise/decay phases (see also 
Montagni et aI7||2006| for a similar behavior in S5 0716+714). 

7 http://users.utu.fi/kani/lm/PG_1424+24Q.html 


Following Danforth et al. ( 2013| >, we modeled these episodes 
with a sum of exponentials after having converted the light- 
curves to flux densities. The rationale is based on the idea that 
the derived time-scales, r, can give constraints on the size of 
the emitting regions. In addition, the time-scales of the decay 
phases, if the emission is due to synchrotron radiation, can allow 
us to derive inferences about the cooling times of the accelerated 
electrons and, in turn, the magnetic fields. 
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The adopted empirical functional form (Danforth et al. 2013| ) 
2F\ 


ex p^) +ex p^y 


reference time scale of 

p4 kx S 


5 min we get R < 3x10 


5 x TfiP c 


10 4 x A cm. The rapid variability identified here amounts to 
only a few percent of the total emitted flux from BL Lac. 

Under the hypothesis that the (variable) emission is due to 
synchrotron and Compton processes, the cooling time-scale can 
limit the time-scale of a decay phase as: 


Ti Cool 


3 m e c( 1 + z) 

4 cr T 6u 0 j e 


s. 


(3) 


Fig. 3. BLLac light-curve after subtraction of the host galaxy contri¬ 
bution and correction for Galactic extinction. A few episodes of rapid 
variability are labelled (see Table[IJ and fits based on Eq. (|T| are also 
shown (blue solid linej. 

Table 1. Parameters of rapid flares during our BL Lac monitoring. The 
epochs are relative to 00:00 UT on 2012 September 2. F, is the ampli¬ 
tude of the variability episode, lcr errors are computed with two param¬ 
eters of interest l |Cas h| 1976j l. 


where m e is the electron mass, cr T the Thomson cross-section, 
u' 0 = u' B + M rad = (1 + q)B 2 1 the co-moving energy density 
of the magnetic field (determining the synchrotron cooling rate) 
plus the radiation field (determining the inverse-Compton cool¬ 
ing rate), q = u. / u' H the Compton dominance parameter, typi¬ 
cally of order of unity for BLLacs (Tavecchio et al. 2010) , and 
y e is the characteristic random Lorentz factor of electrons pro¬ 
ducing the emission. 

The peak frequency of the synchrotron emission is at 


0.274&+B' 

'■\\ n — , " ~ Hz, 


(1 + z)m e c 


(4) 


where e is the electron charge. 

Finally, substituting y e in Eq(s).[3]an the co-moving mag¬ 
netic field can be constrained as: 


(5) 


( 1 ) 


where F\ is the flare normalization, r r ; and Td.i are, respectively, 
the flux rise and decay time-scales, and t, is the time of the pulse 
maximum. The inverse of Eq.[T] is used when the light- 

curve shows a decay followed by a rise, and t { corresponds in 
this case to the pulse minimum. 

The dense sampling of our light-curve allowed us to de¬ 
rive four events with well constrained time-scales (Table[T| and 
Fig.[3]). In all cases the time-scales for rise or decay phases are 
approximately in the range 2-15 min, considering the uncertain¬ 
ties. 

Variability time-scales as short as a few minutes have already 
been singled out for BL Lac objects, mainly at high energies (e.g. 


B > [nm e c{ 1 +z)eloP] ll \;^f c 2 ^5- 112 ~ 

~ 4 x 10 7 (1 + z^C'^ G. 

BLLac is an intensively monitored object. Raiteri et al. 
( j2013] > reported on a comprehensive study of its long-term be¬ 
havior, including the epoch of our observations. From that data 
set the position of the synchrotron peak frequency can be in¬ 
ferred to be close to v syn ~ 5 x 10 14 Hz, and therefore, again 
assuming a reference time scale for decay of ~ 5 min, and con¬ 
sidering that the cooling time should be shorter than this, we get 

B' >6x(^r 1/3 G. 

The SED of BL Lac and that of a number of sources of the 
same class were studied in |Tavecchio et al.| ( |20101 > based on ob¬ 
servations carried out in 2008. A single zone model allowed 
the authors to estimate an average magnetic field B ~ 1.5 G, 
a Doppler factor 6 ~ 15, and radius of the emitting region 
R ~ 7 x 10 4 pc. Compared to the results from our analysis, 
based however on observations carried out in 2012, the emitting 
region of BLLac turns out to be, as expected, a small fraction of 
that responsible for the whole emission and the magnetic field is 


Aharonian et al. 2007; Albert et al. 

20071 [Arlen et al. 2013]) or A similar analysis for PKS 1424+240 is not possible due to 

X-rays (e.g. 

Wagner & Witzel 1995 

), where the flux variation is the very low level of variability shown during our observations. 


a large fraction of the total. In the optical the percentage ampli¬ 
tudes of flux variations are typically lower, possibly due to the 
superposition of several emission episodes with largely differ¬ 
ent time-scales (e.g. Chandra et al.|20lT| |Danforth et al.|[2013| 


Sandrinelli et al. 2014[ > originating from different emitting re¬ 

gions. Therefore, strictly speaking, constraints derived by the 
light-curve analysis hold only for a portion of the emitting re¬ 
gion of the order of the ratio of the flux variability to the total 
flux. 

The size of the emitting region can be constrained as: 


R < 


6ct 
1 +z’ 


( 2 ) 


where z is the source redshift, c the speed of light, and 5 is the 
relativistic Doppler factor of the emitting region. Assuming a 


A fit with a constant is indeed perfectly acceptable although dur¬ 
ing the first ~ 30 min of observations the source was slightly 
brighter by ~ 0.01 - 0.02 mag. 

The length of our monitoring does not allow us to de¬ 
rive general conclusions, although the difference in the ob¬ 
served flux variability between the two objects is remarkable. 
PKS 1424+240 is actually at a higher redshift compared to 
BLLac (z, ~ 0.6 vs. z = 0.069). Time dilation will lead to a 
reduction of any intrinsic variability for the former source by a 
factor of about 1.5 with respect to the latter. In addition, based 
on the SEDs shown in |Tavecchio et al.| ( |2010| l and Aleksic et al.| 
(2014), the optical band is at a higher frequency than the syn¬ 
chrotron peak for BL Lac, and at a lower frequency (or close to) 
for PKS 1424+240. As widely discussed in Kirk et al. ( 1998) , 


under the assumption that magnetic fields in the emitting region 
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are constant, flux and spectral variability depend on the observed 
frequency. If electrons with a given energy, corresponding to 
photons at a given frequency, cool more slowly than they are 
accelerated, variability is smoothed out, as it might be the case 
for PKS 1424+240. Variability is expected to be particularly im¬ 
portant close to frequencies emitted by the highest-energy elec¬ 
trons, where both radiative cooling and acceleration have similar 
timescales. Different short-term variability behaviors for sources 
with the synchrotron peak at lower or higher frequencies than the 
observed band were indeed already singled out (IHeidt & Wagnerl 
[1996][19981 [Romero e t al.|2002||Hovatta et al.|2014). 


In the literature it is also customary to look for the total 
flux doubling/halving times (e.g. Sbarrato et al.|20lT] Impiom 


|bato et ak |2011| |Foschini et al.||2013) . The small amplitude of 

the variability we observed does not allow us to derive strong 
constraints, since this would always require large extrapolation. 
However, the shortest time-scales we could detect are of the or¬ 
der of less than four hours for BLLac, consistent with the values 
found in other blazars, which is consistent with the idea that the 
the whole emitting region is much larger than the regions respon¬ 
sible for the rapid variability. 


A variability analysis can be carried out for the polarimetric 
light curves too. The results show variability timescales at the 
same level as the total flux curves, although with larger uncer¬ 
tainties. Rapid time variability on minute to hours time scales 
for the polarized flux was singled out in other blazars, as for in¬ 
stance AO 0235+164 (Hagen-Thorn et akpOOS), S5 0716+714 
dSasada et al.|2008|), CGRaBS J0211 + 1051 flChandra et al.|2012) 
or CTA 102 ( |Itoh et al.|20T3) . Intranight variability for a set of 


radio-quiet and radio-loud AGN was studied by Villforth et al. 
( [2009) . 


3.2. Polarimetry 


Blazar emission is known to be characterized by some degree of 
polarization that is often variable, both in intensity and direction, 
on various time-scales (see |Falomo et al.|2014[ for a recent re¬ 
view about optical observation of BLLacs). Occasionally, some 
degree of correlation or anticorrelation between the total and po¬ 
larized flux is observed (e.g.|Hagen-Thorn et al.|2008] Raiteri et 


|ak1|2012| |Sorcia et al.||2013[ |Gaur et al.||2014] >, while often no 
clear relation is singled out. The complexity of the observed be¬ 
haviors likely implies that, even when a single zone modeling 
can satisfactorily describe the broad-band SEDs, more emission 


components are actually active. It was proposed (e.g. Barres de 
Almeida et al. 2010] Sakimoto et al. 201 3] that a globally weakly 


polarized fraction of the optical flux is generated in a relatively 
stable jet component, while most of the shorter term variability, 
both in total and polarized flux, originates from the development 
and propagation of shocks in the jet. 


BLLac and PKS 1424+240 show rather different behaviors 
in the linear polarimetry as well. The degree of polarization of 
BLLac starts at about 11% and decreases slowly for a few hours 
to about 9%; then, for the remaining three hours of our monitor¬ 
ing, it decreases more quickly to about 6%. The position angle 
increases rather quickly after the first hour, from about 14° to 
23°; then it remains stable for a couple of hours and then in¬ 
creases again to about 30°. Superposed on these general trends 
there is considerable short-term variability above the observa¬ 
tional errors. PKS 1424+240, on the contrary, shows a fairly con¬ 
stant polarization degree at about 4% and a position angle close 
to 127°, with some variability only at the beginning of our moni¬ 
toring. These behaviors are in general agreement with the results 


reported by Andruchow et al. ( |2005| > studying intra-night polar¬ 
ization variability for a set of BLLac objects. 

The PKS 1424+240 jet was likely in a low activity state, al¬ 
though it was not in its historical minimum (see Sect.[3j. This is 
also confirmed by the publicly available information and data at 
other wavelengths, such as high-energy gamma rays provided by 
the Fermi /LAT Collaboration 8 ] and soft X-rays_available from 
the Swift/X RT monitoring progranj^] Aleksic et al. (2014) re¬ 
ported a higher polarization degree, 7 - 9%, in 2011, when the 
source was brighter than during our monitoring. Lower polariza¬ 
tion degrees, 4.4 - 4.9%, were reported by [Mead et~aL| ( [1990) 
in 1988, when the source was instead fainter. The position an¬ 
gle was about 113 - 119°, similar to that observed during our 
monitoring. The latter is also consistent with the direction of the 
jet as measured by VLB A radio observations at 2 cm (Lister et 
al. 2013). The kinematics of the most robust radio component 
showed a position angle of 141° with a velocity vector direc¬ 
tion of 108°, i.e. with a very small offset (33° + 20°, |Lister et| 
lal.|2013). VLBA observation s in the framework of the MOJAVE 
Projecj^i (Lister et al. 2009) showed a decreasing trend in the 
polarization degree from 5%' in 2011 to 2.8% in 2013, with a 
roughly stable position angle (126° - 154°), consistent with our 
results in the optical. In general, looking at historical data, the 
polarization degree of PKS 1424+240 seems to be almost con¬ 
stant (~ 4%) below a given optical flux (likely < 9.0 mJy, based 
on the refereed studies). The optical position angle seems to be 
quite stable and aligned with the kinematic direction of the radio 
jet and with the radio polarization position angle. This behavior 
might suggest some kind of “magnetic switch” (i.e. a threshold 
effect) in the jet activity (e.g.lPunsly & Coroniti| 19901 [Meier et 
|al.|1997|[T999) . 

Neglecting the short-term variability of BLLac, the total 
rotation of the position angle, taking the minimum at approx¬ 
imately ~ -2 hours (see Fig.[Tj) and the value at the end of 
our monitoring, amounts to about 15°, i.e. 2 - 2.5°/hour (45 - 
60°/day). Rapid position angle rotations of this magnitude are 
not unusual for blazars in general, and for BLLac specifically 
(e.g. [Alter et al. 1981] Sillanpaa et al.||1993| |Marscher et al. 


2008). The observation of relatively stable and long-lasting rota¬ 


tional trends (days to months) suggested that the polarized emis¬ 
sion could be generated in a jet with helical magnetic fields or 
crossed by transverse shock waves, or in a rather stable jet with 
an additional linearly rotating component (Raiteri et al. |2013) . 

Our well-sampled monitoring observations allow us to dis¬ 
entangle different behaviors even during the relatively short- 
duration coverage of BLLac (Fig. [4] upper left plot). At the be¬ 
ginning of our monitoring period, we see a rapid flux decrease 
with polarization slowly decreasing. After that, the source enters 
a phase characterized by rapid small-scale variability both in the 
total and polarized flux. Finally, the flux begins to increase regu¬ 
larly by a small amount and the polarization decreases abruptly 
down to the lowest observed level. The relation between polar¬ 
ization and position angle (Fig. [4] upper right plot) shows the al¬ 
ready mentioned rotation of the position angle with the decrease 
of the linear polarization. However, again superposed on this 
general trend there is considerable variability (see also Hagen- 


|Thorn et al.[2008 


for a similar analysis). 

In Raiteri et al!] (|2013) the long-term (years) flux light curve 


was modeled assuming the flux variation to be (mainly) due to 
Doppler factor variations with a nearly constant Lorentz factor. 


8 http://fermisky.blogspot.it/2®14_06_Ql_archive.html 

9 http://www.swift.psu.edu/monitoring/source.php?source=PKS142 

10 http://www.physics.purdue.edu/astro/MOJAVE/sourcepages/1424+ 
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Fig. 4. (upper left) BL Lac (host galaxy subtracted) flux density vs. linear polarization (host galaxy corrected, assuming unpolarized emission, e.g. 
|Covino et al.|2003] >. At least three different regimes are singled out: at early time the flux changes rapidly with a slowly varying and rather high 
polarization (brown, circles), then an intermediate phase with chaotic flux and polarization variations (green, stars), and finally a sharp decrease in 
polarization with almost constant flux (blue, squares). Times in the legend are in hours (see Fig.|T|. (upper right ) BLLac position angle vs. linear 
polarization. Same symbols of the upper left panel. The position angle tends to increase when the linear polarization decreases. The trend becomes 
very clear at the end of our observation. ( bottom ) Flux density vs. Stokes parameters Q and U. Periods with approximately linear dependence 
between polarization and total flux are also singled out. Same symbols as in the upper left panel. 


i.e. due to small line of sight angle variations. We applied the 
same technique for our rapid monitoring. Knowing the view¬ 
ing angle required to model the flux variations, it is then pos¬ 
sible to predict the expected polarization in different scenarios. 
In the case of helical magnetic fields, following Lyutikov et al.j 
(20050 . we can derive a polarized flux fraction at 9-10%, roughly 
in agreement with our observations. However, a detailed agree¬ 


ment, explaining the short-time variability for both the total and 
polarized flux, is not possible. Alternatively, we may consider 
transverse shock wave models ( |Hughes et al.|1985) >, with which 
again rough agreement for the polarization degree is reached, but 
no detailed agreement is possible. A geometric model for the flux 
variation is therefore unable to simultaneously interpret the total 
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flux and polarization behavior at the time resolution discussed 
here. 

As already introduced in Sect. 3.1 a possible interpretation 
of both total and polarized flux curves can be derived if it is as¬ 
sumed that the observed emission is due to a constant (within 
the time-scale of our monitoring) component with some degree 
of polarization and one (or many, Brindle et al. 1985) rapidly 
varying emission component(s) with different polarization de¬ 
gree and position angle (see also Sasada et al. 2008; Sakimoto et 


al.|2013| >. The idea is rather simple; using the first three Stokes 


parameters the observed polarization can be described as: 



Abs — 

Aonst~l~ 

Amt 


s = ■ 

SobsAbs — 

(2const Aonst 

0 var Aar 

(6) 


Aobs Abs — 

U const Aonst 

^AarAar 



where the suffixes "obs", "const" and "var" refer to the observed 
(total), constant and variable quantities. The redundancy in Eq.[6] 
can be reduced following various possible assumptions, often 
depending on the availability of multi-wavelength datasets or 
long-term monitoring (see, e.g.|Holmes et al. 1984||Quian|1993[ 


Brin dle||1996 Barres de Almeida et al.||2010| l. Hagen-Thorn et 


al.|( 20021 ) assumed, based on their long-term polarimetric mon¬ 


itoring, that the stable component in BLLac could be character¬ 
ized by P ~ 9.2% and 0 ~ 24°. As discussed in Hagen-Thorn 
|& Marchenko| ( |1999| > and |Hagen-Thom et aT] ( |2008| l, if a lin- 
ear relation between polarized and total flux is singled out, this 
can allow one to estimate the polarization degree and position 
angle of the variable component. A linear relation between the 
Stokes parameters and the total flux implies that polarization de¬ 


gree and position angle are essentially constant (Hagen-Thorn & 
|Marchenko| 1999[ > and their values can be derived as the slopes 
of the linear relations. At the beginning of our monitoring we 
can identify a sufficiently long and well defined linear relation 
between the Stokes parameters Q and U and the total flux (see 
Fig[4] bottom panel). As already mentioned, we find consider¬ 
able variability superposed on the linear trend. Neglecting the 
shorter term variability, we can roughly estimate P rdr ~ 22% 
and 8 ~ 34°. The constant component turns out to be remarkably 
consistent with the one identified by Hagen-Thorn et al. ( |2002[ l 
at a flux level ~ 9.5 mJy. 

At any rate, the increasing complexity singled out by long- 
and short-term monitoring requires new theoretical frameworks 
for a proper interpretation. Marscher ( 2014[ >, for instance, pro¬ 
posed a scenario in which a turbulent plasma is flowing at rela¬ 
tivistic speeds crossing a standing conical shock. In this model, 
total and polarized flux variations are due to a continuous noise 
process rather than by specific events such as explosive energy 
injection at the basis of the jet. The superposition of ordered and 
turbulent magnetic field components can easily explain random 
fluctuations superposed on a more stable pattern, without requir¬ 
ing a direct correlation between total and polarized flux. As dis¬ 
cussed in |Marscher| ( ]2014| , simulations based on this scenario 
can also give continuous and relatively smooth position angle 
changes as observed during our monitoring of BLLac. 

[Zhang et al.| ( |20l4j ) presented a detailed analysis of a shock- 
in-jet model assuming a helical magnetic field throughout the 
jet. They considered several different mechanisms for which a 
relativistic shock propagating through a jet may produce a syn¬ 
chrotron (and high-energy) flare. They find that, together with a 
correlation between synchrotron and synchrotron self-Compton 
flaring, substantial variability in the polarization degree and po¬ 
sition angle, including large position angle rotations, are possi¬ 
ble. This scenario assumes a cylindrical geometry for the emit¬ 


Table 2. Parameters of interest for the model based on the scenario ex¬ 
tensively discussed in Zhang et al. (2014) and Zhang et al. ( 2015j). An¬ 
gles are in the co-moving frame. 


Parameter 

Value 

Bulk Lorentz factor 

15 

Length of the disturbance (L) 

3.8 x 10 14 cm 

Radius of the disturbance (A) 

4.0 x 10 15 cm 

Orientation of the line of sight 

90° 

Helical magnetic field strength 

2.5 G 

Helical pitch angle 

47° 

Electron density 

4.5 x 10 2 cm -3 


ting region moving along the jet, which is pervaded by a helical 
magnetic field and a turbulent component. On its trajectory, it 
encounters a flat stationary disturbance, which could be a shock. 
This shock region does not occupy the entire layer of the emit¬ 
ting region, but only a part of it. In the comoving frame of the 
emitting region, this shock will travel through the emitting re¬ 
gion, and temporarily enhance the particle acceleration, resulting 
in a small flare. After the shock moves out, the particle distribu¬ 
tion will revert to its initial condition due to cooling and escape. 

The 3DPol (3D Multi-Zone Synchrotron Polarization) 
code presented in |Zhang et al.| ( |2014|) and MCF P (Monte 
Carlo/Fokker-Planck) code presented in |Chen et al. (201 ^re¬ 
alizes the above model. As elaborated in Zhang et al. ( j2015| , 
since the shock is relatively weak and localized, the enhanced 
acceleration will lead to a small time-symmetric perturbation in 
the polarization signatures. Some of the key parameters for the 
model are reported in Table 2] and the fits to the polarization de¬ 
gree and position angle light curves are shown in Fig. [5] Near 
the end of the observation, the polarization degree experienced 
a sudden drop, while the position angle continued to evolve in 
a time-symmetric pattern. Therefore an increase in the turbulent 
contribution is necessary although, due to the lack of a multi¬ 
wavelength SED, it cannot be well constrained. The total flux, 
given the very low variability amplitude observed during our 
observations, was set at a constant level of about 12mJy (see 
Fig. [3]). Nevertheless, rapid polarimetry clearly reveals its diag¬ 
nostic power, showing need of inhomogeneity and turbulence in 
the emitting region. 


4. Conclusions 

In this work we are presenting results from rapid time- 
resolved observations in the r band for two blazars: BLLac and 
PKS 1424+240. The observations were carried out at the 3.6 m 
TNG and allowed us to measure linear polarimetry and photom¬ 
etry almost continuously for several hours for both sources. In 
practice, long-term monitoring observations of relatively bright 
blazars can only be achieved with dedicated small-size tele¬ 
scopes; however the richness of information obtainable with a 
rather large facility as the TNG allows us to study regimes which 
were in the past only partially explored. 

BLLac and PKS 1424+240 show remarkably different vari¬ 
ability levels, with the former characterized by intense variabil¬ 
ity at a few percent level, while the latter was almost constant 
for the whole duration of our observations. The shortest well 
constrained variability time scales for BLLac are as short as a 
few minutes, allowing us to derive constraints on the physical 
size and magnetic fields of the source regions responsible for the 
variability. 
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Time (hours, 0 = 2012, Sept 2nd 0:0 UT) 


Fig. 5. Fit of the Zhang et al. ( 2014, 2015 ) model scenario (green solid 
line) to the BL Lac linear polarization (host galaxy corrected, assuming 
unpolarized emission, e.g. |Covino et al.|2003| l and position angle data. 
The general behavior is fairly well described by the model, although 
for the polarization the addition of weakly constrained turbulence in the 
magnetic field is required. 

The variability time-scales for the polarization of BLLac 
are compatible with those derived for the total flux, while 
PKS 1424+240 shows an almost constant behavior also in the 
polarization. The position angle of BLLacertae rotates quasi- 
monotonically during our observations, and an analysis of the 
total vs. polarized flux shows that different regimes are present 
even at the shortest time-scales. 

Different recipes to interpret the polarimetric observations 
are considered. In general, with the simplest geometrical models, 
only the average level of polarization can be correctly predicted. 
More complex scenarios involving some turbulence in the mag¬ 
netic fields are required, and promising results are derived by 
a numerical analysis carried out following the framework de¬ 
scribed in Zha ng et~aT7| ( |2014| |20151 >, which requires some sym¬ 
metry in the emitting region, as shown by the time-symmetric 
position angle profile. The time-asymmetric polarization profile, 
and its decrease during the second part of the event, which is 
accompanied by a few small flares, can be described by adding 
some turbulent magnetic field structure to the model. 
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